Bessel processes, the integral of geometric 
Brownian motion, and Asian options 

Peter Carr and Michael Schroder 



Abstract. This paper is motivated by questions about averages of stochastic pro- 
cesses which originate in mathematical finance, originally in connection with valuing 
the so-called Asian options. Starting with [Y], these questions about exponential 
functionals of Brownian motion have been studied in terms of Bessel processes us- 
ing the Hartman- Watson theory of [Y80] . Consequences of this approach for valuing 
Asian options proper have been spelled out in [GY] whose Laplace transform results 
were in fact regarded as a noted advance. Unfortunately, a number of difficulties 
with the key results of this last paper have surfaced which are now addressed in 
this paper. One of them in particular is of a principal nature and originates with 
the Hartman- Watson approach itself: this approach is in general applicable without 
modifications only if it does not involve Bessel processes of negative indices. The 
main mathematical contribution of this paper is the developement of three principal 
ways to overcome these restrictions, in particular by merging stochastics and complex 
analysis in what seems a novel way, and the discussion of their consequences for the 
valuation of Asian options proper. 

Keywords. Asian options, integral of geometric Brownian motion, Bessel processes, 
Laplace transform, complex analytic methods in stochastics. 

1. Introduction: This paper addresses questions about exponential functionals of 
Brownian motion and the integral of geometric Brownian motion in particular. These 
questions reduce to the study of the quadratic variation processes A^'^^ of geometric 
Brownian motion which for any real drift v are explicitly given by the integrals over 
time 



with B a standard Brownian motion. These processes have both a surprisingly rich the- 
ory and manifold applications ranging from the physics of random media to mathemati- 
cal finance and insurance. In fact, the insurance motivated study of certain perpetuities 
in [Du90] seems to have initiated this line of research. Here, the above integrals over the 
whole time axis are considered and shown to be distributed as the reciprocals of certain 
gamma variables. Drawing on his probabilistic interpretation of the Hartman- Watson 
identities in [Y80], Yor was able to extend this work and to determine the law of the 
processes A^'^^ in [Y]. This approach, using the Laplace transform and based on Bessel 
processes, has been found by Yor to open many surprising vistas on the processes A^'^^ 
and their applications, see in particular [YOl]. The interest in these processes AM IS 
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partially due to their importance in mathematical finance, in particular for understand- 
ing the so-called Asian options. 

Asian options are widely used financial derivatives. As discussed in Part I of the paper, 
these options provide in general nonlinear payoffs on the arithmetic average of the 
price of an underlying asset. A common objective in their valuation is to derive an 
explicit expression for a certain functional of A^'^\ The pursuit of this objective has 
evolved over the last fifteen years as an interplay between theoretical and computational 
perspectives, see [RS] or [DuOO] for instance. Yor's work in [Y] clarified the structure of 
the Black-Scholes prices of Asian options by expressing them as certain triple integrals. 
In contrast to this result, however, it was the Laplace transform approach of [GY], also 
based on Hartman- Watson theory, which had fax reaching consequences for the way 
Asian option valuation is seen today, and on which we focus here from Part II onwards. 
This is essentially because a very explicit expression in terms of Kummcr's confluent 
hyper geometric function resulted from this approach for what has been regarded as 
the Laplace transform of the value of such options. Prom a numerical point of view, 
this expression has proved to be amenable to computation by numerical inversion, 
sec [FMW] as a recent example. And it seems fair to say that this has lead financial 
mathematics to a new interest in developing and applying these techniques. 

Unfortunately, some difficulties with this Laplace transform approach to valuing Asian 
options have emerged. First, it turned out that the Laplace transforms computed in 
[GY] are not those of the Asian option's value. This appeared to empty the relevance 
of this result for the finance application proper. Luckily, it turned out that there is 
a reduction of the original problem of valuing Asian options to the one considered in 
[GY]. All of this is discussed in Part II of the paper. 

A difficulty of a more serious and more principal nature, however, originates with the 
Hartman- Watson approach on which the Laplace transform computations of [GY] are 
based. Its idea is to analyze A('^) using Bessel processes of indices u, and its applicability 
has limits if u is negative because of the pathologies which Bessel processes of such 
negative indices develop. Together with background material about Bessel processes, 
we have thus given in Part III a new exposition of the analysis in [GY], which expli- 
citly takes care of the nonnegativity restriction on the index u. Hereby, we have been 
encouraged by the kind support of Yor, and we have tried to incorporate his kind 
tutorials and suggestions. In terms of financial mathematics, our extension to negative 
ly in fact extends the analysis of Asian options from the zero dividend situation originally 
considered in [GY] to one of general real dividend yields, and thus to one with general 
real risk-neutral drifts. In this setting, the condition that u is not negative transcribes 
into the postulate that risk-neutral drift is not less than half the squared volatility. 
Unfortunately, this lower bound on the drift restricts the financial applicability of the 
results. In fact, the greater is the volatility, the greater is the range of parameters in 
which the nonnegativity condition on u is violated and in which the approach does not 
give the Laplace transforms. Unfortunately, it is precisely due to high volatility of the 
underlying asset that Asian options are used in the first place. 
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Thus, the third contribution of our paper is to remove these restrictions by showing 
the existence of the desired Laplace transforms. In fact, we develop three ways for 
removing these restrictions. More precisely, we develop three princpal ways of coping 
with the difficulties caused by Bessel processes of negative indices in the Hartman- 
Watson approach. This can be seen as the main mathematical contribution of the 
paper. The approaches of Parts IV and VI arc based on an analysis of Bessel processes 
in the spirit of Yor. Our first approach extends the [GY] analysis to the missing cases 
using deeper properties of Bessel processes of negative indices. Our second approach 
gives a new uniform proof using zero index Bessel processes only. Our third approach 
of Part V extends that of [CS] and merges stochastic with complex analytic techniques. 
This appears to be a rather novel and promising line of attack as we are able to largely 
dispense with Bessel processes and focus instead on Brownian motion. 

Apart from all this, our extension of the Laplace transform approach of [GY] has made 
possible advances in valuing Asian options, some of which are sketched in Part VII. 
Hence, it may be fair to say that the Laplace transform approach of [GY] has proved 
to be a rich source for new results and insights in both finance and mathematics. 

I Prologue 

2. Black- Scholes modelling: The results to be discussed originate from the so- 
called risk-neutral approach to the valuation of contingent claims. General equilibrium 
treatments for this and other notions developed for the analysis of financial markets and 
instruments are in [D88], [D96], [KS, Chapters 1-4], for instance. This analysis is based 
on models of security markets, and this section aims to recall the most fundamental of 
these, the Black-Scholes model of security markets. 

In fact, we need only that particular case of the Black-Scholes model where there are 
only two securities, and the understanding is that these arc traded on markets where 
their prices are determined by equating demand and supply. First, there is a riskless 
security, a bond, whose price (3 grows at the continuously compounding positive interest 
rate r, i.e., for which we have we have /?< = expert) at any time t G [0, oo). Then, there 
is a risky security. The fundamental idea is that all uncertainties affecting its price S 
yield a certain probability space. In fact, consider for this a complete probability space 
equipped with the standard filtration of a standard Brownian motion on the time set 
[0, oo). Giving expression to the fact that S comes as an equilibrium price, we have the 
risk neutral measure Q on this filtered space, a probability measure equivalent to the 
given one. And with B any standard Q-Brownian motion, the exact modelling then is 
that S is the strong solution of the following stochastic differential equation: 

dSt = w-Sfdt + a-Sf dBt, t G [0, oo) , 

or equivalently using Ito calculus, St = Sq ex.p{{w—a'^/2)t+aBt), for all t in [0, oo). The 
positive constant a is the volatility of S. The specific form of the otherwise arbitrary 
constant zu depends on the nature of the security modelled which could be a stock, a 
currency, a commodity etc. For example, if S is the price of a stock paying a dividend 
continuously so as to have constant dividend yield 5, then we have zu = r—S. 
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3. Asian options and their equilibrium pricing: In the Black-Scholes framework 
of §2, fix any time to and consider the process J given for any time t by: 

ft 



J{t) = / Su du , 

J to 



The arithmetic- average Asian option written at time to with maturity T and strike 
price K is then the stochastic process on the closed time interval from to to T paying 

^(^^ KV = ma./0,^-i. 



T-to J { T-to 

at time T and paying nothing at all other times. As such it is a contingent claim on 
the time interval from to to T with payoff (J(T)/(T — to) — K)^ . 

It is one of the fundamental insights that in the the equilibrium framework of the 
Black-Scholes model, any such contingent claim on a risky security has an equilibrium 
price which is equal to the discounted expectation of its payoff with respect to the risk 
neutral measure conditional on today's information, see [D88, §22K, (47)], [D96, §8A] 
or [MR, Corollary 5.1.1], for instance. Applying this arbitrage pricing principle, the 
price Ct of the Asian option at any time t between to and T is given by the discounted 
Q-expectation conditional on the information ^ available at time t: 



a = e-'^^-'^E'^ 



T-to 



■S^t 



However, following [GY, §3.2], we do not focus on this price, but instead we normalize 
the valuation problem as follows. On factoring out the reciprocal of the length T— to 
of the time period, we split the integral J{T) into two integrals, one of which is deter- 
ministic by time t and the other of which is random. We then couple the deterministic 
integral with the new strike. For the random integral, we restart the Brownian motion 
driving the underlying at time t, and then using the scaling property of Brownian mo- 
tion, we change time so as to normalize its coefficient in the new time scale to two. The 
precise result is the factorization: 

T-to 

which reduces the general valuation problem to computing 

the normalized time-t price of the Asian option. Herein, A^'^^ is Yor's process 



Ct='-,^-^-C(^Hh,q), 



and the normalized parameters are as follows: 

2tU , , (7^ 

cr^ 4 



h=—{T-t), q = kh+q*, 
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where 

K ( /"* 

A: = — , q = q{t) = —^\K ■ (t-to) - Sudu 

To interpret these quantities, u is the normalized adjusted interest rate, h is the nor- 
malized time to maturity, which is non-negative, and q is the normxilized strike price. 
On a conceptual level, notice that valuing any Asian option in this way is reduced to 
computing a single function C'^'^\ and that a similar notion of normalized hedging of 
Asian options can be developed along these lines. On a structural level, moreover notice 
how q becomes affinc linear in the time variable h with coefficients k and q* depending 
only on quantities known at time t. 

4. A first reduction of tlie normalized valuation: There is now a dichotomy in 
computing the normalized time-i price 

C^'^\h,q)=EQ[{A^;:^-qf] 

of the Asian option according to the normalized strike price q being positive or not. 
Indeed, if q is non-positive, Asian options lose their option feature. Computing the 
values C^^\h,q) is straightforward: 



Lemma: If q <0, we have 

C^''\h,q)=EQ 

with 



A^"^ -a 

g2/i(j/+l) _ ]_ 



for any real v, and it is thus sufficient to compute C^'^\h, q) if q > 0. 

This can be proved on applying Fubini's theorem, and the last expectation is seen to be 
analytic in v with its value at v = —1 equal to h. It should be mentioned that formulas 
for all moments of A^'^^ have been derived at various instances over the last fifty years, 
see for example [YOl, §2.4.1, (4.d'), p.33 and Postscript #3b), p.54]. 

5. Yor's integral representation for Asian option values: A closed form for the 
normalized time-t prices C^*^) of Asian options can be obtained as a consequence of 
Yor's triple integral representation [Y, (6.e), p. 528]. Recall the latter is based on Yor's 
characterization of the law of A^'^^ in [Y, (6.c), p. 527] and so is eventually based on the 
Hartman- Watson theory of [Y80]. Furnishing a measure for the difficulty of computing 
normalized prices, the precise form of Yor's closed form is as follows 

Theorem: For any reals h, q > and u, we have 

Ci-){h,q)=c,,h X- e ■ (--q) ■^^yih)dydx. 

Jo Jo ' 
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Herein the function for any positive real number a, is given by the following integral 

-acosh(u>) . , , . / TT 



e 27i g 



sinh('u;) • siny—wj dw , 



,2h 



Jo 

for any h > 0, and we abbreviate 

1 

TTyZTT'^h 

While Yor's formula seems to require v to be bigger than at least minus one, it is valid 
for all real v. This is proved in [SL, §6]. 

One purpose of closed form expressions is to provide means for actually computing 
option prices. While Yor's results above are the key to many insights into the math- 
ematical structure of A^'^\ Yor's formula of the Theorem has a number of structural 
difficulties in this regard. First, it involves three integrations with seemingly no further 
structure or further possibilities for simplification. Thus, methods for explicitly com- 
puting it will necessarily be rather complex. However, we have noticed an apparently 
even bigger obstacle to computability. For instance taking to = and zu = r equal to 
5%, we compute the factors Ci,^h as follows: 





a = 20% 


a = 30% 


a = 40% 


T = 1 year 


2 . 627 X 10213 


2 . 265 X 10'^'' 


4.816 X 10^2 


T = 6 months 


7 . 686 X 10^^27 


5.717 X 10^^^ 


2 . 583 X 10^°^^ 



Table 1. c„,h as function of t and <t. 

As the examples of §8 or §24 will show, normalized prices of Asian options are not 
too big. Yor's formula thus expresses them as the product of a big number times a 
triple integral. The latter so has to be very small and must be computed with very 
high accuracies to get reasonably accurate results. The Laplace transform approach 
developed in [GY] explicitly for the purpose of valuing Asian options was seen to offer 
a way out of these difficulties. 



II Laplace transform results 

6. The Laplace transform in option valuation: Working with continuous func- 
tions on the non-negative real line of at most exponential growth, the Laplace transform 
S'{f) of any such function / is defined by 

/•oo 

S'{f){z)= / e— /(x)dx, 
Jo 

for any complex number z in a half-plane contained sufficiently deep within the complex 
right half-plane. 

The connection of this notion with option valuation in general and Asian option valua- 
tion in particular is as follows. Fix any option type, like §3's European style Asian call 
option on a certain stock with price S. It will depend on a number of parameters, like 
strike price and maturity date. At a fixed point in time t, consider the family which 
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consists of all options on the market with all such parameters fixed except maturity 
dates M. In the Asian option example thus consider the Asian options of all maturities 
available at time t which have the same strike price. In this way regard maturity date 
M as a real variable ranging from t to infinity. The value of the option thus becomes a 
function of M. 

However, it is normalized prices C^'^^ we have to consider for the Asian option. The 
normalizations of §3 in fact turn C*-^^ into a function of, in particular, normalized time 
to maturity. As a function of maturity date M normalized time to maturity is explicitly 
given by h{M) = {cr'^/4:) ■ (M — t). With M from t to infinity, h{M) thus ranges from 
to infinity. The normalized price of the Asian option so becomes a function on the 
non-negative real line. Call it /ao ^or the sake of emphasis. Recalling from §3 how 
the normalized strike price q depends in an affine linear way on normalized time to 
maturity, we then have more precisely 

fAo{x)=E'^[{A^:^-{kx+q*)y' 

for any non- negative real x. 

At this point we should signal a difficulty with [GY] which we have noticed in Fall 1999. 
It is not the functions /ao this last paper is working with, but the functions fGY,a given 
for any positive real a by 

fGYA^) = EQ[{Ai''^-a) 

for any non- negative real x. Notice that the function /ao has a non-constant strike 
price X kx + q* whereas any function fGY,a has the constant strike price a. So the 
value function of the Asian option /ao is different from any hmction fGY,a which thus 
is the value function of a non-Asian option. Because of the injectivity of the Laplace 
transform on continuous functions, the Laplace transform 5='(/^o) is different from any 
Laplace transform .S"{fGY,a) too. The problem at this point so is one of relevance of the 
mathematics for the finance application proper: is there a way of relating the valuation 
of Asian options to the valuation of the non- Asian options? 

7. Valuing Asian options using families of non-Asian options: The basic idea 
is as follows. Try to reconstruct the normalized time-t price 

Ci-){h,q) = E'^Ua^^'^ -{kx + q*))^] 

i x=h 

of the Asian option from a family of auxiliary functions whose single members are 
unrelated to the problem of valuing the Asian option but amenable to the [GY] ana- 
lysis. Given that the time dependency of the normalized strike price poses the problems, 
simply force this strike price to be constant. Thus arrive, for any positive real a, at §6's 
functions fGY,a on the positive real line recalled to be given by 

/Gy,a(x) = S'3[(AM-a)+] 

for any positive real x. These are the functions considered in [GY]. As we remarked 
in §6 they are the values of certain non-Asian options, and taken individually, they 
cannot be used to value the original Asian option. However, our finding is that as a 



8 



P. Carr and M. Schroder 



whole they allow one to recover the normalized time-i price. With the concepts of §3 
this key reduction is made precise in our following 

Lemma: If q = kh+q* is positive, computing the normalized time-t price C^'^\h,q) of 
the Asian option reduces to computing all fGY,a with a > 0. More precisely, C^^\h,q) 
is obtained by choosing the function fGY,kh+q* ^-^'^ evaluating it at h. 

A moment's reflection will convince the reader that this is true by construction, and 
meanwhile, the particular case of the Lemma we explained to Yor in May 2000 can 
be found in [YOl, pp. 95-96]. Notice that our construction works in the more general 
situation where functions of the form h{y) = f(y, (p{y)) with a known map (p have to 
be computed: compute the functions / and then intersect with the graph of if to get h. 
And to stress the main result again, as one consequence of this technique, we have at 
this point reduced valuing Asian options to valuing the family of all non- Asian options. 

8. Laplace transforms of the non- Asian option values fGY,a- The key reduction 
of the preceding §7 shows the way to apply the [GY] Laplace transform to value Asian 
options. Moreover, adopting the notation of §3, consider the family of all functions 
fGY,a with a > that send any a; > to 

/Gy,a(^) = ^^[(4^)-a)+]. 

Recall from §7 that its single members are unrelated to valuing the Asian option, 
but that as a whole, the family allows reconstruction of the normalized time-t price 
C(^\h,q) of the Asian option. As a first step in actually computing the fGY,a, try to 
compute their Laplace transform FGY,a given by 

/•oo 

FgyA^) = / e-^VGy,a(^) dx = S'{fGY,a){z). 
Jo 

Here the complex number z is to be taken in a half-plane sufficiently deep within the 
right complex half-plane such that the integrals are finite. The function so obtained 
is analytic. The precise conditions under which these integrals are finite is part of our 
description of these generalized Geman- Yor Laplace transforms in the following 

Theorem: // the normalized strike price q is positive, the integrals FGY,a o,re finite for 
any complex number z with Re (z) > max{0, 2{i'+l)}, and we have 

D^{a,z) 
^^^'"^"^ = z{z-2{u+l)) 

where on choosing the principal branch of the logarithm 

1 

\dx . 



Herein is the modified Bessel function with complex order fi, as discussed in [L, 
Chapter 5]. Generalizing [GY, (3.9), p. 363], the integral Dy{a^z) can be expressed 
using the confluent hypergeometric function $ discussed in [L, Chapter 9] as follows 
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Corollary: For any complex z with Re (z) > max{0, 2(i/+l)}, we have 
on setting ijl{z) = V^z+i^. 

We again stress that these Laplace transforms are not those of the Asian option price, 
but rather the Laplace transforms of the prices of auxiliary options. To obtain Asian 
option prices, we have to invert these Laplace transforms and then proceed using 
§7 Lemma; formally speaking: 

C^''\h,q) = S^^-'{FGY,,){h). 

In full mathematical generality, analytic inversion has been achieved in [SA], and we 
come back to this in §23. Numerical inversions have also been accomplished. For exam- 
ple [FMW] computed the following seven cases as reproduced in [DuOO, Table 7.1] 



Case- 


r 


a 


T 


■S'o 




2C(-) 


1 


2% 


10% 


1 


2.0 


3 


0.056 


2 


18% 


30% 


1 


2.0 


3 


0.219 


3 


1.25% 


25% 


2 


2.0 


-0.6 


0.172 


4 


5% 


50% 


1 


1.9 


-0.6 


0.194 


5 


5% 


50% 


1 


2.0 


-0.6 


0.247 


6 


5% 


50% 


1 


2.1 


-0.6 


0.307 


7 


5% 


50% 


2 


2.0 


-0.6 


0.352 



Table 2. Prices 2C<"> for k=2.o and 
to=t=o using numerical Laplace inversion. 



Obtaining our two results proceeds in two steps. In a first probabilistic step, the argu- 
ments of [GY] are adapted to compute the modified Geman-Yor transforms FGY,a- We 
give in Part III a new exposition of the argument incorporating the tutorials and kind 
suggestions of Yor. The key idea is to factorize the geometric Brownian motion of the 
underlying over a Bessel process of index ly. Pertinent notions are discussed in §9. This 
makes time stochastic in such a way that Yor's process now takes the double role 
of both a stochastic clock and a control variable for the Asian option. At first sight, 
this appears to complicate the original valuation problem. However, this double role of 
A^'^^ is especially suited to the Laplace transform. Indeed, in contrast to the situation 
for the Asian option, the strike price a of the non-Asian option with value function 
fGY,a is independent of time. This makes it possible to reduce the computation of the 
Laplace transform Fcv^a to the following problem: obtain explicit expressions for the 
Bessel semigroup of index v and for a certain conditional expectation involving first 
passage times of Yor's process A^'^^ . 

However, such results are available for both concepts only if the index v is non-negative, 
an assumption which is explicit in [GY, §2]. This non-negativity condition, however, 
translates into the condition that the risk-neutral drift is not less than half the squared 
volatility. Unfortunately, this places restrictions on the financial applicability of the 



10 



p. Carr and M. Schroder 



result. For example, if volatility is 30%, the arguments of [GY, §3] are not valid if the 
difference between the riskfree rate and the dividend yield is less than 4.5%, and then 
we do not have the Laplace transforms of the non- Asian options either. Worse yet, the 
greater is the volatility, the greater is this range of parameters in which we do not have 
the Laplace transforms. Unfortunately, it is precisely due to high volatility that Asian 
options are used in the first place. 

We have been a bit disconcerted by these findings. Luckily, however, we found that 
in particular those results of Table 2 where u is negative were reproduced in [DuOO] 
using an alternative approach. And corroborating Andre Weil's [We, p. 457] dictum that 
'theorems are proved by those who believe in them', we are now able to discuss in the 
sequel three different ways for stablishing the Theorem and its Corollary for arbitrary 
real risk neutral drifts u. 

The first of these, as discussed in Part IV, is inspired by [Y] , and we think Yor could 
have given this argument had he been aware of the financial motivation for letting the 
parameter u be negative. In fact, the key idea is to try to bypass the difficulties of 
Bessel processes of negative index i/ by Girsanov transforming to the simpler Bessel 
processes of index zero.This theme is developed also in the third approach discussed in 
Part VI. Here, the idea is not to Girsanov transform Bessel processes derived from the 
geometric Brownian motion driving the underlying. Instead, Girsanov transform this 
geometric Brownian motion itself; an idea we distilled from [Y80]. The effect of this 
is that zero drift Bessel processes enter right from the beginnning, and the result is a 
uniform argument based on these most natural Bessel processes. 

In comparison to these two approaches, our second approach discussed in Part V seems 
somewhat novel. The idea is to combine stochastic methods with complex analytic ones. 
The net effect here is that with an input of some standard result from the latter area 
like the identity theorem, it is not Bessel processes which now have to be dealt with 
but Brownian motion. For this it is moreover not required to work on the process level, 
as it is STifficient to work on the level of expectations. So this second approach seems 
to be an example of a rather promising methodology for solving problems, which is to 
systematically enhance stochastic techniques with complex analytic ones. 



9. Preliminaries on Bessel processes: As a preliminary to estabhshing the Laplace 
transforms of §8, this section collects a number of pertinent facts from the theory of 
Bessel processes. This theory is patterned after the example of the Bessel processes of 
integer dimension 5 > 2, which are defined by taking the Euclidean distance from the 
origin of a Brownian motion in dimension S. Applying Ito's Lemma, their infinitesimal 
generator .jy is seen to be given by 



III Laplace transforms if z/ > 



2u+l 
2x 
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for any function / in C|(R>o). This notion makes sense for any real number 6, and the 
real-valued diffusion associated to using the Volkonskii construction, see for instance 
[K, Theorem 4.3.3, p. 91], is the Bessel processes i?H on [0, oo) with index v = ((5/2) -1. 
While Bessel processes of non-negative indices v stay positive if started with a positive 
value at time zero, Bessel processes of negative indices v develop some pathologies. As 
explained in [RY, XI §1], in this case they hit zero. If —1 < i/ < 0, they are thereupon 
instantantanously reflected and never become negative. For v = —\ they continue at 
zero. 

This matters for the second way of defining Bessel processes of arbitrary dimension b. 
In fact, the focus here is on squares of Bessel processes. Applying Ito's Lemma, they 
are the continuous strong solutions of the stochastic differential equation 



[RY, XI, §1]. These stochastic differential equations make sense for any real number 
(5 and have a unique continuous strong solution also if 8 is smaller than two. The so 
obtained processes are studied in [YG6, §3]. For non- negative indices z^, they coincide 
with the squares of the corresponding Bessel processes of index v. They develop some 
pathologies for negative indices v. In this case, they hit zero if started with a positive 
value at time zero, and if < — 1, they even continue negative. Notice that in such 
situations their square roots are purely imaginary, and so cannot coincide with any 
of the Bessel processes constructed above. However, these two ways of extending the 
notion of Bessel processes do coincide for processes started at time zero at a positive 
value up to the first time zero is hit. This essentially is the reason behind the following 
Lamperti identity, which may nevertheless be surprising 

Lemma: For the index v any real number, we have the factorization: 



for any t > 0, where — e^'^^'^'^'^^^dw is Yor's process. 

For u > a proof is given in [Y92, §2] while the general case is now contained as 
exercise XI (1.28), p. 452 in the third edition of [RY]. We are indebted to Yor for this 
and for kindly supplying us with the following argument. 

The idea is to apply the Ito rule to the square Zt of Yt = exp{ut+Bt) to get 



Time change the process using the inverse function T{t) = 'mf{u\J^ Z^j dw > t} to Yor's 



dpt = 25dt + 2 





process to get 




w ■ 
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To interpret the stochastic integral in this sum, apply the basic time change formalism 
for stochastic processes as in [0, §8.5] to obtain 

Jo 

where Wt is defined as the stochastic integral Wt = f^^^'' \/Z^ dBy, and is a Brownian 
motion. Using the inverse function theorem of calculus, the derivative of r is equal 
to the reciprocal of the derivative with respect to time of Yor's process A^'^^ at time 
w. Hence it is equal to the reciprocal of Zy,. On substitution we so identify the time 
changed process Z as a continuous solution to the stochastic differential equation for 
the square of the Bessel process of index u: 

=2(i/+l)t + 2 ( ^/z^dW^. 
Jo 

Using the uniqueness of the solution of these stochastic differential equations, the time- 
changed process Z is the square of the Bessel process of index u. Reversing the time 
change, this translates into 

To establish the identity of the Lemma, we have to take square roots. This is not 
a problem if is non-negative since then the Bessel process takes non-negative values 
only. It does pose a problem if v is negative. In this case, however, recall that the Bessel 
process starts at time zero with a positive value. Since it is continuous by hypothesis, 
it will stay positive until it first hits zero at time t* > 0. Since the process A^'^'^ starts 
at zero at time zero, there is a latest point in time t**, infinity admitted, such that A^'^') 
is smaller than t* at all points in time t smaller than t**. Thus we have the required 
identity at least for all points in time t smaller than t**. Now Yt is never zero. Since 
the processes on both sides of the identity are continuous in time, t** must be infinity, 
and the proof is complete. 

10. Computing Laplace transforms if z/ > 0: This section is the first step in the 
proof of the integral representation of §8 Theorem for the Laplace transform 

/•oo 

FGY,aiz)= / e~^''fGY,a{x)dx = S'{fGY,a){z), 

Jo 

where, with the concepts of §3 and §6, we have fGY,a{x) = E^[{A'^^'> —a)'^], for any pos- 
itive real numbers a and x. We now explain why one needs to restrict the probabilistic 
arguments of [GY] and apply them mutatis mutandis in order to arrive at the following 

Lemma: The assertions of ^8 Theorem are valid if v = 2a~'^w—l > 0. 

We are very indebted to Yor for correspondence and discussions about this result, and 
are very grateful for his kind support. In the sequel, we want to indicate the key steps 
of the proof following [GY], while trying to incorporate his suggestions. Hopefully, no 
pitfalls have remained undetected. 
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The basic idea is to make time stochastic using the Lamperti identity 

for all positive real numbers w as has been discussed in the preceding section. Here, i?*^*^^ 
is the Bessel process of index v with i2(^)(0) = 1. On applying this Lamperti identity, 
A^"") has the double role of both control variable and stochastic clock. That the "strike 
price" a is independent of time now becomes essential. It makes possible to transcribe 
the condition on the control variable to be bigger than a as the inverse time change 

T^^a = inf{^t I A^u"* > a} 
for the stochastic clock. This is the key idea for obtaining the representation 

J2{u+l)lw-T.,a]+ _ I 



fGY,aH = 



2 



for all w > 0. Indeed, fix any positive real number x, and consider the process 
at X on the set of all events where Ti,^a takes values less than or equal to x. Break the 
integral defining A^'^'^{x) at Tn^a- The first summand then is A^'^^ at time Tj/^o and so 
is equal to a. In the second summand, restart the Brownian motion in the exponent 
of the integrand at Tj/^a shifting the variable of integration accordingly. The second 
integral then is the product of exp(2 • {B{Ti,^a)+i^Tu,a)) times Ai") at x— Tj^^a, by abuse 
of language after having applied the Strong Markov property. This last process is such 
that it is independent of the information at time Ti,^a- Unravelling the definition of r^^a, 
the first above factor so is the square of the Bessel process at time a. Now taking 
the expectation conditional on the information at Tj/^a, we thus get: 



EQ\{AP-ay 



A 



On substitution for the expectation of A^'^\w) from §4 Lemma or using [Y, §4], the 
required expression for fGY,a follows. 

At first sight this appears to complicate the problem. However, it is just what is espe- 
cially suited to the Laplace transform FcY^a of fGY,a now given by: 



FgyM) 



1 



Mn2 



dw . 



Still, for computing this integral one wants to interchange the Laplace integral with 
the expectation . If z is real, it seems best to follow Yor's proposal for justifying 
this. Indeed, with the integrand of the double integral in question positive and mea- 
surable, apply Tonelli's theorem now but justify only in a later step that any of the 
resulting integrals are finite. The case of a general argument z is reduced to this case 
by considering the absolute value of the integrand, and the result is the identity 



FgyM 



z{z-2{u+l)) 
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of measurable functions for any complex number z with Re {z) > 2(i/+l). The idea for 
identifying the expectation in the numerator as Di,{a,z) then is to condition on the 
Bessel process to obtain 

/•oo _ , 

Du{a,z)= / x^E^ e-^^"''' = X puA'^,x) dx , 
.70 L J ' 

where pi,,a is the time-a semigroup density of the Bessel process of index ly starting at 1 
at time zero. Following [GY, p. 362] we make this integral explicit by making the single 
factors of its integrand explicit. For this, work with the results recalled in [GY, §2]. 
With respect to the conditional expectation factor, under the hypothesis > 0, Yor 
has computed it at positive real arguments z in [Y80, Theoreme 4.7, p. 80] (see also 
[GY, Lemma 2.1 and Proposition 2.6]). Using analytic continuation, the validity of his 
result can be seen to extend to the arguments z in the right half-plane required in the 
present situation. This then gives for the conditional expectation factor in D,^{a, z) the 
following expression as a quotient of I-Bessel functions: 

eQ L---" it;W(a) = J = ^-^^^^ {-) . 

Explicit expressions for the Bessel semigroups p^^a have been known for > — 1 for 
some time, see [Y80, (4.3), p. 78] or [GY, Proposition 2.2]. The density pi,^a{^,w) of the 
time-a Bessel semigroup with index v and starting point 1 is 



Pi/,a(l, W) = e '^a . _ . 

a V a / 



Nothing seemed to have been known about such densities if z/ < —1 before [YG6]. 
However, the results proved there for // < — 1 still need to be handled with care as it 
will be explained in §14. The upshot is that, in accordance with [GY, §2], the above 
decomposition of Dy{a,z) seems to give explicit results without further qualifications 
only if if z/ > 0. Then, however, we have the required result 

D,{a,z) = -^J^ e-^x^+'l^^[-)dx. 

There is a further technical point to be taken care of herein: choose the principal branch 
of the logarithm to define the square root on the complex plane with the non-positive 
real line deleted. 

To complete the Tonelli argument proposed to us by Yor and to complete the proof, 
we have to establish the finiteness of this integral for any fixed complex number z with 
Re (2:) > 2(z/ + l). This is a consequence of §11 Proposition's convergence analysis of 
these integrals, and granting this result, the proof of the Lemma is complete. 

11. Integrability analysis: In terms of the concepts of §3, this section studies finite- 
ness of the integrals 

1 

D,{a,z) = ^J^ e-^x''+^I^^[-)dx 
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for any real a > in terms of their complex parameters v and z. The precise result is 
the following 

Proposition: Let e > 6e any real number. If |Im(z/)| < e, the integrals Di,[a,z) are 
finite for any complex z with real part Re (z) > 2e^ . 

The Proposition depends on the following result about the complex square root asso- 
ciated to the principal branch of the complex logarithm. 

Lenrnia: Let e > be any real number. For any complex v with |Im(^')| < e we then 
have 

Re {V2z + z/2 ) > |Re(z/)| 
for any complex z with Re (z) > . 

A proof of the Lemma based on a close analysis of the square root can be found in 
[CS, §10]. To prove the Proposition, finiteness of D^(a,z) under the conditions of the 
Proposition follows by combining the above square root lemma with the asymptotic 
behaviour of the Bessel function factor of its integrand near the origin and towards 
infinity. Indeed, setting 



= V2z+u^ , 



from [L, §5.11] recall that is a continuous function on the positive real line in par- 
ticular whose asymptotic behaviour for large real arguments is 

Hence the factor exp(— x2/(2a)) dominates the asymptotic behaviour of the integrand 
of Di,{a,z) with x to infinity, whence its integrability away from the origin. On the 
other hand, from [L, §5.7] we have for real arguments near zero 

Thus, if the real part of ji+v+A is positive, or cquivalcntly, if we have 

Re(/x) > -(Re(z/)+4), 

no integrability problems arise for x near the origin. Under the conditions of the Propo- 
sition, on the other hand, the above Lemma gives 

Re(|Lt) > |Re(i/)|. 



Since this last inequality implies the former, the proof of the Proposition is complete. 
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IV Laplace transforms in the general case: 
using Girsanov transforms of Bessel processes 

12. Further preliminaries on Bessel processes: Our first way of extending the 
results of §10 to negative indices v requires further preliminaries on Bessel processes 
from [Y, §2]. Recall that Bessel processes where the index v is any real number are the 
real-valued diffusions whose infinitesimal generators .sif are given by 

^/(x) = lr(x) + ^/'(:r), 

for any function / in C^(R>o). To describe the law P^l,■u, C(R>o,R>o) of i?^'^) if this 
process starts at the non-negative real u, let p be the canonical process on C(R>o, R>o); 
recall it operates as evaluation map: /9a(/) = f{p)- If ^ is the canonical filtration with 
•/^ a equal to the sigma algebra generated by the pg with s < a, we then have the 
mutual absolute continuity relation 

Lemma: // the Bessel process of any real index v is started at any non-negative real 
u, its law is realated to that of the zero index Bessel process started at u as follows 

^^^'^\ .^an{a<To} -\uj 2 Jo pV '""l ' 

where Tq is the first passage time of p to zero. 

This is proved as an application of Girsanov's theorem by exchanging drifts in the 
stochastic differential equation of §9, and has the following 

Corollary: For any complex z with positive real part, and any non-negative real r, 
where the expectation is taken with respect to the law Po,u- 



The Corollary is proved in two steps. If z is any non-negative real, Tq = oo, and we 
have the explicit expressions for the densities Pn,a{u,r) of the of the Bessel semigroups 
already encountered in §10 and related to the law via P^^u(a, dr) = P/^,a(n, r) dr: 

P„a{u,r) = exp (- Yaiu'+r^))l.{^), 

for any non-negative reals u, a > and r. And so the Corollary follows on taking 
expectations in the absolute continuity relation of the Lemma. Observing that both 
sides of the identity to be proved are analytic functions in z on the right half plane, 
the general case then follows by analytic continuation as a second step. 
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13. First proof of the Laplace transform using Bessel processes: Resuming 
the discussion of §10, we are now able to complete the proof of §8 Theorem in the way 

it might have been envisaged by Yor: based on a careful analysis of Bessel processes. 
Recall that we still need explicitly compute for negative normalized risk-neutral drifts 
u, the risk-neutral expectations 



where a > and Re (z) is positive and sufficiently big in particular. From the time 
change part of the argument in §9, recall that T,y^a as the inverse time change of the 
process A('^) at time a is given by 

ds 

' v,a 



Again using §9's Lamperti relation 



(i?ro 



which is valid for w > 0, the point now is that a is smaller than the first passage time 
to zero To of the Bessel process R('^\ Put differently, this Bessel process which starts 
in 1 at time zero still is positive at time a. Using the absolute continuity relation of 
§12 Lemma, 



R 



pie 



^rds 

Ps 



r ds-i 



Conditioning on the index-0 Bessel process thus gives 



Jo 



El 



ds 



Pa 



Ps 



/5''"^^P0,a(l,/') dp, 



With the explicit form of the semigroup recalled in §12, applying §12 Corollary gives 

.2 



as desired. An application of §11 Proposition then shows finiteness of this integral if 

zq = Re (z) is positive and bigger than 2(z/ + l), and the first proof of §8 Theorem is 
complete. Note that meanwhile the argument of[Y01, pp. 97-99] seems to corroborate 
our statement that this proof is very much in the spirit of Yor. 

14. Vista on the use of Bessel processes: The strategy of the preceding argument 
is to bypass the difliculties brought about by Bessel processes with negative indices by 
Girsanov transforming to Bessel processes of index zero. One could ask about a direct 
attack in the spirit of §10. This would require possession of explicit expressions for 
both the densities of the respective Bessel semigroups and the pertinent conditional 
expectations of exp(— zTi^^a). The recent work of Going- Jaeschke and Yor in [YG6] now 
in particular provides certain analytic expressions for the densities. However, it is not 
Bessel processes as considered in §12 which are studied there, but rather processes 
obtained as strong solutions to the stochastic differential equations of squared Bessel 
processes, as mentioned in §9. Recall that for negative indices, the latter processes 
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become negative and their square roots, which should give the Bessel processes, then 
are purely imaginary. Still, the two notions of Bessel processes thus obtained coincide 
on their respective positive range. There, we have Bessel semigroup densities based on 
those derived in [YGo, §3 Proposition 2, p. 21]. For indices v < —1 and time-0 starting 
values y > 0, they are given by 

Pu,t{x,y) = h{x,y,d,t)e J — e dw 

defining 

h(xv5t) = - i^yy^fS-s 

and with 5 = 2(1 + < and n = l — v. However, these densities are in terms of 
new classes of special functions. The clarification of their relations to those for Bessel 
processes of non-negative indices is but one of the problems that require further study. 



V Laplace transforms in the general case: 
combining stochastics and complex analysis 

15. Remarks on general philosphy: The discussion up to now has focused in 

particular on extension at the process level. The actual valuation problem, however, is 
not at the process level but at the expectation level. From this point of view, §10 has 
identified two function / and g in the variable u and has proved them to be equal for u 
non-negative. One would like to have this equality also for negative and thus extend 
the validity of the identity f = g from the non-negative real line to the whole real line. 
Such situations quite commonly appear in problems in analysis and are addressed there 
using analytic continuation. However, functions become amenable to complex analytic 
methods only on open subsets of the complex plane. Thus the identity theorem of 
complex analysis asserts that two functions on a connected open subset of the complex 
plane are equal if they are analytic and agree on a convergent sequence there only. So 
it is in fact no longer possible to stick to real numbers only. As a subset of the complex 
plane they arc closed with an empty interior. At this stage, however, nothing is known 
about the functions of §10 if v is outside the non-negative real line. In particular, 
it is not known if they exist at all, and this needs to be established togother with 
their analyticity properties. Of these two functions, Dy{a,z) is already given as an 
explicit analytic expression, while the other function is not. In fact, it is not explicit 
at all, as it is defined as the Laplace transform of a certain expectation. The question 
to be tackled here is then how to get explicit analyticity properties from such non- 
explicit stochastic concepts. In the present situation, we attack this in two stages. 
First, establish analyticity properties of the expectation. Then, in a second step, study 
how these are preserved on taking Laplace tranforms. As it turns out, using standard 
results from complex analysis, we find that it is not Bessel precesses which enter into 
the analysis but simply Brownian motion. All of this may be regarded as an instance 
for why enhancing stochastics by complex analytic methods seems quite an interesting 
and promising line of thought. 
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16. First step - analyticity of the function Di,{a, z): As a first step in tlie ana- 
lytic continuation argument, this section studies the analytic properties of the general- 
ized first Wcbcr integral Dy{a, z) of §8 Theorem. On choosing the square root associated 
to the principal branch of the logarithm recall 

is finite for any positive real number a and for any complex numbers z and u such that 
the real part of i/ + 4 + (2z + i^^)^/^ is positive as a consequence of §11 Proposition's 
integrability analysis. Using the confluent hyper geometric function $ discussed in [L, 
§9.9], the precise analyticity result to be proved is the following 



Proposition: Let a be any positive real and e be any non-negative real. For any complex 
number v with |Im < e we then have 

if z is any complex with Re (z) > setting fi = s/^z+u^. 



Corollary: Let a be any positive real and e be any non-negative real. For any complex 
number z with Re (z) > 2e^, sending v to D^{a,z) then defines an analytic map on the 
set of all complex numbers u with |Im {y)\ < e. 

Both results are based on §11 Proposition which gives finiteness of D^{a, z) under their 
conditions on u, a, and z. Combining this with the analyticity properties of $ discussed 
in [L, §9.9] and those of the gamma function, the Corollary follows from the Proposition. 
The proof of the Proposition then reduces to exlicitly computing D^{a, z). For this we 
modify the quite typical discussion in [W, §13.3, pp.393f] of Hankel's generalization 
of Weber's first integral. The idea is to expand the modified Bessel function in the 
integrand of 




into its series of [L, §5.7] and integrate term by term. Using [L, §9.9] this is justified by 
the absolute convergence of the series for the confluent hypergeometric series which is 
to result, and we get 

(2a)'^;fr^r(/x+l+n) n! 
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Changing variables y = (2a) ^x'^, compute any n-th integral as 

Extracting the series of the pertinent confluent hypergeometric function we thus get 



: — . cb I — '■ — -;- (/ -J- 1 • — \ . On 

r(/x+i) 

Multiplying this expression with exp(— (2a)'~^)/a, the Proposition follows. 



17. Second step - Analyticity of the functions fGY,a' Establishing analyticity 
results in v about the Laplace transform of the expectation defining non- Asian option 
prices combines insights from stochastics with insights of an analytic nature. We estab- 
lish this result in two steps. As a first step, in this section any of the auxiliary functions 
fGY,a of §7 is considered as function in the variable u 

L{x,v)=EQ[{AP-af] 

for any fixed positive real numbers a and x. Using §10 Lemma, we know it is defined 

for non-negative real numbers v. However, this has been achieved in a very indirect 
way only: for these values of v the Laplace transforms of the corresponding functions 
fGY,a have been shown to be finite. Now more is true indeed 

Lemma: For any x > the function v ^ L{x, v) extends to a function on the complex 
plane which is analytic at each point, and for which we have the majorization 



For the proof of the Lemma now set f{w) = {w — a)~^. Applying Girsanov's theorem 
such that Wx = vx+B^ becomes a standard Brownian motion, and dropping reference 
to this new measure, we get 



L{x,v)=E\f{Af)e^ 



X,,2 



For establishing the analyticity statement of the Lemma, it is thus sufficient to show 
that the expectation factor is analytic in any complex number v. This is true by defi- 
nition if we have the convergent series 



E 



m=0 



for all v. For this it is sufficient to show that the series is absolutely convergent for all 
u. Using the Cauchy-Schwarz inequality this is implied by the convergence of 



ml 



E 



lit: V 

m=0 

for all u. Herein E[f'^{A^x^)\ is majorized by the second moment of Yor's zero drift 
process at x, and so is finite from [Y, §4]. Since the factors ^'[W^"*] are majorized 



E 



W2" 
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by 7r~^/^ • (2a;)™' • ml for all m > 0, convergence follows using the ratio test. Actually, 
we have so established yet another upper bound to the price of the Asian option. 

To establish the majorization of the Lemma, taking absolute values inside the expec- 
tation in the above Girsanov representation of L{x, u) gives: 



\L{x,u)\<E\f{Af) 



e 



The absolute value of the exponential factors are the exponentials of the real parts of 
the respective arguments. Majorizing the function in Yor's zero drift process by this 
process itself, we get 



,-fRe2(,.) 



Reversing the Girsanov transformation then completes the proof of the Lemma. 

18. Third step - analyticity of the transforms F^y „: While the preceding sec- 
tion has studied the expectations 

L{x,v)=E'i[{A^^)-af] 

for any fixed positive real numbers a and x as function in the complex variable v only, 
this section moreover treats x as a variable. If u is any non-negative real number, we 
have from §10 Lemma that the integrals of the Laplace transforms 



/•oo 
JO 



F{v){z) = / e-^''L{x,u)dx 



are finite if Re (2;) > 2(z^ + l). In this section we give an independent proof of the 
following more general statement 

Proposition: For any complex number z with a positive real part, the map sending u 
to F{i>){z) is analytic in all complex numbers v with Re (z) > ^Im^(z/)+2(Re (z^) + l). 

The proof of the Proposition is based on the following 

Lemma: For any complex number v, the Laplace transform F{v){z) is finite for any 
complex number z with Re {z) > max{0, |lm ^(i/) + 2(Re (i/) + l)}. 

Proof of the Lemma: As first step in proving the Lemma, we establish for any complex 
number v the majorization 



Jo 



for any complex number z in the sense of measurable functions. Indeed, majorize the 
absolute value of F{i'){z) by taking the absolute value inside the defining integral. The 
absolute value of the exponential factor then is equal to exp(— Re {z)x). Majorizing the 
absolute value of L(x,z^) using §17 Lemma, the estimate follows. 
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Setting vq = Re (i^), now let Re (z) be positive and such that = Re {z)—lm (i^)^/2 is 
bigger than 2{uo + l). We compute the Laplace transform of the right hand side of the 
above inequality using §4 Lemma. If uq is different from minus one, 

roo -I 

Jo ?o(?o - ^(t'o + ijj 

using that is bigger than 2(i/o+l) to compute the improper integrals. It converges to 
^q"^ with going to minus one. Thus it is seen to coincide with the Laplace transform 
for the case uq = —1. The Laplace transforms are finite if is positive and bigger than 
2(fo+l), and then F{i/){z) is finite a forteriori. The proof of the Lemma is complete. 

Proof of the Proposition: For proving the Proposition, fix any complex z with a posi- 
tive real part and then choose any complex z/q satisfying the rcsTilting inequality of the 
Proposition. In fact, the validity of this inequality then extends to all in a compact 
neighbourhood V of vq. As a consequence of §17 Lemma, fz{x, u) = exp(— v) is 
for any x > analytic in i/ on F in particular. Since V is compact, the argument of the 
Lemma moreover shows that its absolute value is majorized by an intcgrable function 
g on the real line. The function u F(y){z) so is continuous on V being obtained by 
integration of fz over the positive real line. To show it is analytic on the interior of V 
we want to apply Morera's theorem, see [R, 10.17, p. 208] and have to show 

/ F{v){z)du = ^ 
JdA 

for any triangle A contained in the interior of V. Indeed, since we have shown fz to be 
integrable, applying Fubini's theorem gives 

/ F{u){z)di'= / / fz{x,u) du dx . 
JdA Jo JdA 

Recalling that i' i— > f^[x,u) is analytic from §17 Lemma, the inner integral herein is 
zero by Cauchy's theorem for A, see [R, 10.13, p. 205]. Thus the whole double integral 
is zero, as was to be shown. This completes the proof of the Proposition 

19. Final step - second proof of the Laplace transform using analytic exten- 
sion: It remains to pull things together and show how the results of this Part V give 
a second proof of the two results of §8 by using analytic extension. 

First notice that §8 Corollary is implied by §8 Theorem using the computation of 
Dy{a,z) in §16 Lemma. Thus we are reduced to prove §8 Theorem. 

The proof of §8 Theorem is by analytic continuation using the identity theorem, see 
[R, Corollary to 10.18, p. 209]. As a consequence of §10 Lemma it remains to establish 
the identity in §8 Theorem for negative indices u only. The existence of the Laplace 
transform on all complex numbers z with positive real part required in the Theorem 
then has essentially been proved in §18 Proposition. 

For establishing the crucial Laplace transform identity of §8 Theorem, thus first let 
z be any complex number with real part Re (z) > 4 and choose < e < 1. Using 
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§16 Corollary, the generalized first Weber integral 

1 

is analytic in f on the e-thickened real line which consists of all complex numbers 
u with |Im {u)\ < e. Picture as the band of height 2e symmetric with respect to the 
real axis. If we choose e so small that Re {z) > 2(2+2e), we claim to have analyticity 
of the Laplace transform 

f oo 



poo 

F{u){z)= / e-'''L{x,i^)dx 
Jo 



10 

as a function in v on the £-thickened half- line Bg,, that subset of A^ which consists of 



/ 





















Figure 1. The e-thickened half-line Be. 

all complex numbers v with |Im {u)\ < e and Re (ly) < e. Indeed, if Re (i^) < e, we have 
2(2-1- 2£) > 2£-F2(Re(i/) + 2). If |Im(z/)| < e, we have 2e > lm'^{iy)/2 since e < 1. Any 
u in so satisfies the inequality Re {z) > Im^(z/)/2+2(Re (i/)+2) of §18 Proposition, 
and the claim follows. 

So the functions F{i/){z) and 

Du{a, z) 



Dl{a,z) 



z{z-2{u+l)) 

are analytic as functions in v on the e-tliickcncd half-line Bg,. It is now a consequence of 
§10 Lemma that we have F{v){z) = D*{a,z) for all > such that 2(z/+l) < Re(z). 
With Re (z) > 4 this so holds a forteriori for all > such that 2(i^ + l) < 4, i.e., 
for all non-negative real numbers u smaller than 1. With e smaller than 1 we so have 
F(i'){z) = D*{a, z) for any u in the subinterval (0, e) of S^. With Bg, open and connected 
the identity theorem literally applies to give that F{i')(z) = D*{a,z) for all v in B^. 
This identity then holds a forteriori for all real numbers u in B^, i.e., for all v < e. 
This gives §8 Theorem for any z with Re {z) > 4. 

To lift this last restriction on Re (z), notice that §18 Lemma implies F{u) for any fixed 

real v to be analytic on the half-plane {Re (z) > 2{i' + l)}; this is seen by using a 
Morera type argument as for proving §18 Proposition. On the other hand, applying 
§16 Corollary, D*{a, z) as a function of z is analytic on the intersection of this last half- 
plane with the right-half plane. The validity of the identity F{u){z) = D*{a^ z) thus can 
be analytically continued from complex numbers z with Re {z) > 4 to complex numbers 
z with Re {z) positive and bigger than 2(i/+l), and the second proof of §8 Theorem is 
complete. 
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VI Laplace transforms in the general case: 
a uniform proof 

20. Remarks on general philosophy: The question may arise if there is not a 
uniform way for computing the Laplace transfoms in §8 Theorem not requiring a two 
step procedure. The Lamperti factorization 

of §9 Lemma may in fact offer a clue. The idea of the two arguments discussed up to 
now is to focus on the Bessel process side of this identity. It is the problems with Bcssel 
processes that transcribe into their two-step approach. A possible remedy so would be 
to try to bring in Bessel processes not at the earliest stage but as late as possible. The 
Lamperti factorization suggests to focus on the geometric Brownian motion of its left 
hand side. While Yor has already made extensive use of this device for instance in [Y92, 
§6], we explain in the sequel what seems to be a most natural adaption of the Girsanov 
technique to the Lamperti factorization approach for explicitly computing the Laplace 
transforms FGY,a with a > 0. 

21. Thh"d proof of the theorem - a uniform argument: This section sketches a 

uniform way for computing the Laplace transforms Fcy^a with a any positive real, and 
thus provides a third proof of §8 Theorem. Recall these transforms are defined by 

/•oo 

FGY,a{^)= e-^''fGY,a{x)dx 

Jo 

for any complex z with Re (z) sufficiently big, and where 

/Gy,a(x) = i?«[/(AM(x))]. 

setting f{x) = (x — a)+. The key idea is to apply Girsanov's theorem such that Wx = 
vx + Bx becomes a Brownian motion, and suppressing reference to this new measure, 
transcribe fGY,a in terms of Yor's zero drift process A(^) as follows 



/Gy,a(x) = i^«[/(4°V"^ 

With To, a = inf{n I ^(f^) (u) > 0} the inverse time change to A^^^ at time zero, the 
computations of §10 thus prove fGY,a{x) = on the set of all events where x < To,a- 
And on the set of all events where x > ro,a we now have the representation 



(T}(0)\2 _ ^(0) 



f{A(o)) 

which applies to express fGY,a{x) as the following iterated expecation 



/Gy,„(x) = E 



4(0) 



Applying Laplace transforms to both sides of this identity then gives 



FgyA^) = F 



.JO 



To, a 
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In the conditional expectation restart the Brownian motion at time ro,a- To the resulting 
additional v-th power of exp(WT-o,a) apply the corresponding Lamperti factorization and 
interpret it as being equal to the i^-th power of the zero index Bessel process at time 
a. Collecting powers of this last process, 



(0)..+2 (.+ ^) 



TO, 



r 

/ 

Jo 



To, a 



The Strong Markov condition which we have used here also makes the whole inte- 
grand in the conditional expectation independent of time-ro,a information. Reversing 
the Girsanov transformation in this inner expectation, 



FgyA^) = E 



{R^^^y+'e<'^^y°-E\ r e-i'+^^^AP dx 

I Jo 



This puts us into the situation of §10. Partially reversing the Tonelli argument, the 
inner expectation is equal to the first moment of the drift-i^ process A('^) at time x 
recalled in §4 Lemma. Computing the Laplace transform, we thus arrive at 



1 



1 



E 



(0)xl/+2 -(z+i^)TO,a 



2(z/+l) Vl-2(i/+l) 

if Re {z) is bigger than at least max{0, 2(i^+l)}. In this way we are reduced to compute 
the expectation factor. This, however, proceeds as in §10 by conditioning on the cor- 
responding Bessel process. The whole point is that this time conditioning is not on a 
Bessel process of an arbitrary index but on a Bessel processes of index zero. Observing 
how u enters now via the exponent of the index zero Bessel process and via a shifting 
factor for the time change ro,a, we so arrive at 



E 



1 

'2a 



dx , 



a JO 

for any z with sufficiently big positive real part. Using §11 Proposition, this integral is 
finite if Re {z) is positive and bigger than 2{i'+l). So this third proof of §8 Theorem is 
complete. 



VII Epilogue 

22. Consequences for Asian options - Hermite functions: Having persevered 
to this point, the reader may wonder about the nature and the quality of the impli- 
cations of the mathematics developed up to now. Indeed, going back to the starting 
point of the journey, it seems that the Laplace transform approach makes possible 
significant improvements in understanding §3's normalized prices C^''\h,q) of Asian 
options themselves. These improvements are on a structural level and as such make 
possible advances on computing C^^^ as a consequence. This is essentially by being 
able to establish new links of Asian option valuation with a well-studied class of special 
functions: the Hermite functions to be reviewed in this section. 

Following [L, §§10.2ff], to which we refer for details, Hermite functions Hfj_ are analytic 
on the complex plane as functions of both their variable z and their degree jj,. If the 
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real part Re (//) of n is bigger than —1, they have the integral representation 
Hi^{z) = J x^cos \ 2zx—-iJ,Trj dx . 

So they specialize to the fi-th Hermite polynomials if /x is any non-negative integer, 
whence Hq = 1, Hi{z) = 2z, H2{z) = 4z'^-2, H^{z) = 8z^-122;, for example. If the real 
part of n is negative, however, Hermite functions change their character. Then they 
have the integral representation: 



1 



and so specialize via {2/ ^)H_i{z) = exp(2;^)Erfc (z) to the complementary error 
function Erfc recalled to be given by 



2 f°° t2 

Erfc {z) = ^ e'^ . 



For any complex Hermite functions can be expressed in terms of the Kummer con- 
fluent hypergeometric function $ by 

^^-r(V2^ /_M 1 2'^-T{-l/2) fl-, 3 ,^ 

' r((l-/x)/2) I 2'2'^>'+'^ r(-/x/2) I 2 '2'^ J 

for any complex z. Prom this representation one can moreover see how Hermite functions 
arc connected with the parabolic cylinder functions -D^ and with the Kummer confluent 
hypergeometric function of the second kind ^. 

23. Consequences for Asian options - new integral representations: Hermite 
functions as recalled in the previous section appear naturally in the closed form solu- 
tion for §3's normalized price C^^\h, q) of the Asian option we have developed in [SA]. 
If the normalized strike price q is positive, it expresses this value as the sum of inte- 
gral representations which have a product structure. They are obtained by integrating 
products of Hermite functions with weighted error functions as follows: 

Theorem: // q is positive, the normalized price C^^'Xh^q) of the A sian option is given 
by the following difference 

where the are three-term sums 

= Qrig^(p^) + Chyp,^(p5) + Chyp -^(p^) 

whose single summands are integrals that depend on parameters > 0, but which as a 
whole are independent of these. 

In terms of §3's concepts, c is given by 

, , r(i.+4)-(2g)5(-+2) 

c = c(z/, q) = J — 

27r • (i^+l) • 6^9 
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recalling v = 2w ja'^—V. With p any non-negative real, the trigonometric terms Ctrig€(p) 
with ^ equal to or z/-|-2 are the integrals 

over the real parts of products of Hermite functions times certain functions Eij{h), and 
the hyperbolic terms Chyp,^(p) with ^ equal to ±1^ or ±(1^+2) are the integrals 

over the imaginary parts of such products. Herein E^{h) are the weighted complemen- 
tary error functions for any complex w given by 

E^ihKw) = e^Wc (-^ + . 
Remark: If p equals zero, the trigonometric terms specialize to 

Comparing with the formula of §5, the above formula, is given as a sum of single 
integrals whose integrands have a structural interpretation as products of two functions. 
It identifies the higher transcendental functions occuring as factors in these products, 
and shows how they are given by or built up from Hermite functions. 

24. Epilogue: On a technical level, the differences just noted between §5's and §23's 
formula can be regarded as consequences of the different mathematical approaches for 
proving the valuation formula. In fact, §5's originates from Yor's direct attack on the 
law of the integral of geometric Brownian motion. In contrast, §23's is eventually based 

on the indirect enveloping construction of §7. In fact, to obtain Asian option prices, 
first analytically invert §8 Theorem's Laplace transforms FGY,a, then proceed using the 
key reduction of §7 Lemma; formally speaking: C^'^\h,q) = 2f~^{FGY,q){h). 

However, there are not only structural differences between §5's and §23's formula. Fi- 
nally deriving benchmarks for the normalized prices C^^\h,q) of Asian options seems 
to be one of the main practical application of such formulas. As we have already men- 
tioned, Yor's formula seems impracticable for this purpose in particular because of the 
gigantic size of the numbers it involves. For instance, consider valuing Asian options 
with annual interest rates r equal to nine percent, with maturities of one year, with 
K = Sq, with to = 0, and with a volatility of cr = 30%. These values require coping with 
numbers of order 10^'^''. Our formula improves on this aspect as well. In the above situ- 
ation, for instance, the four hyperbolic terms have orders 10^ and the two trigonometric 
terms have orders 10~^ . Sharpening results of [RS] , it thus became possible to derive in 
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[SH, Chapter 5] the following first time benchmark vahies for normalized Asian option 
prices 



a 


maximal error 




20% 


4.9727 X 10"^'' 


. 00074155998788343 


30% 


4 . 9687 X IQ-^*^ 


. 00217354504625037 


40% 


4.9157 X lO-^*^ 


. 00478100328341654 


50% 


4.9461 X 10-^*^ 


. 00890942045213227 



Table 3. Normalized prices C^''\h,q) of the Asian option for t=i. 

All of this may be seen as a consequence of the Laplace transform approach to valuing 
non- Asian options initiated in [GY]. With hindsight, this paper appears to be a rich 
source for new results in both mathematics and finance. 
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